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The stabilized jellium model (SJM) provides us a method to calculate the volume changes of 
different simple metals as a function of the spin polarization, £, of the delocalized valence electrons. 
Our calculations show that for bulk metals, the equilibrium Wigner-Seitz (WS) radius, r s (C), is 
always an increasing function of the polarization i.e., the volume of a bulk metal always increases as 
C increases, and the rate of increasing is higher for higher electron density metals. Using the SJM 
along with the local spin density approximation, we have also calculated the equilibrium WS radius, 
f s (N, C) , of spherical jellium clusters, at which the pressure on the cluster with given numbers of total 
electrons, N, and their spin configuration C vanishes. Our calculations for Cs, Na, and Al clusters 
show that f 3 (iV, C) as a function of C behaves differently depending on whether N corresponds to a 
closed-shell or an open-shell cluster. For a closed-shell cluster, it is an increasing function of C over 
the whole range < C < 1, whereas in open-shell clusters it has a decreasing behavior over the range 
< C < Co, where Co is a polarization that the cluster has a configuration consistent with Hund's 
first rule. The results show that for all neutral clusters with ground state spin configuration, Co, the 
inequality r s (N, Co) < ^s(O) always holds (self- compression) but, at some polarization Ci > Co, the 
inequality changes the direction (self-expansion). However, the inequality f s (N, C) < ^s(C) always 
holds and the equality is achieved in the limit N — ► oo. 



I. INTRODUCTION 



Jellium -model ( JM) with spherical geometry is the simplest model used in theoretical study of simple metal 
clusters. EIlI In the spherical JM, the ions are smeared into a uniform positive charge background sphere of den- 
sity n = 3/47rrg and radius R — (zA r ) 1 / 3 r s where z and N are the valence of the atom and the number of constituent 
atoms of the cluster, respectively. In the simple JM calculations the value of r s is taken to be the bulk value of the 
Wigner-Seitz (WS) radius of the metal. In some calculations. involving metal surfaces, one may-.use diffuse JM to take 
account of the ion relaxations at the surface of that metaljj However, it is a well-known factum that the simple JM 
yields negative surface energies at high electron densities (r s < 2), and negative bulk moduli for r s as 6. Also, using 
the simple JM, it is not possible to find a realistic size and energetics for metal clusters. These drawbacks have roots 
in the mechanical instabilityQtJ of the simple JM system. That is, the bulk jellium system is stable only for r s — 4.18. 
In 1990, Perdew et al. have introduced the stabilized jelliumEl model (SJM) by adding two corrections to the simple 
JM energy. The first correction subtracts the spurious self-energy of each WS cell from the JM energy, and the second 
correction adds to the energy the effect of difference in the potential an electron sees from the discrete pscudoions and 
from the jellium background. The second correction introduces the core radius parameter, r c , of the pseudopotential 
which can be adjusted in such a way that the stability of the bulk jellium system could be achieved at any observed r s 
value for the valence electrons. Using the value for the core radius ( which stabilizes the unpolarized bulk system) 
in the SJM energy functional of a spin-polarized bulk system or a cluster, it is possible to find the equilibrium r s 
value for any given spin polarization Q. 

In this paper we have calculated the changes in the equilibrium r s values due to polarization for bulk Cs, K, Na, 
Li, Ga, Al, and for Cs, Na, Al clusters using the SJM in the framework of local spin density approximation (LSDA). 
Our calculations show that for the bulk system the volume is always an increasing function of the polarization and 
the increasing rate is higher for higher electron density metals. Also, we have solved the self-consistent Kohn-Sham 
(KS) equations!! for Cs, Na, and Al clusters of different sizes in the spherical SJM and obtained the equilibrium sizes 
of the clusters for each spin configuration consistent with Pauli's exclusion principle. The results show that in a 
closed-shell cluster with N total electrons and polarization £, the equilibrium WS radius f s (N,() is an increasing 
function of C over the whole range < C < 1 whereas, for an open-shell cluster (except for the two nearest neighbors 
of a closed-shell cluster), it has a decreasing behavior over < £ < Co, an d an increasing behavior over Co < C < 1- 
Here, Co corresponds to an electronic spin configuration for which Hund's first rule is satisfied. It has been shown 
previouslyES that this property also holds in the case of the simple JM with spherical geometry. Also, the equilibrium 
r s values for a neutral cluster are found to be always smaller than that of the bulk metal with the same value of the 
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polarization This effect is due to the surface tension which is appreciable for small clusters. On the other hand, 
for an N electron neutral cluster there exists a polarization value Ci j beyond which. thp equilibrium r s value of the 
cluster exceeds the bulk r s value of unpolarized metal. This is called self-expansionllilil2l and for two different neutral 
clusters with the same number of electrons N, the £1 value of the higher electron density metal is smaller than that 
of the lower density metal. 

The structure of this paper is as follows. In section [n] we first show the method by which one can calculate the 
changes in r s values due to spin polarization for different bulk metals. We then explain how to apply the method to 
metal clusters in order to calculate the equilibrium r s values for different spin configurations. In section III we present 
the results of our calculations and finally, the work is concluded in section |y|. 



II. CALCULATIONAL SCHEME 



In the context of the SJM, the average energy per valence electron in the bulk with density parameter r s and 
polarization £ is given byiLj 



£{r s ,C r c) = t s (r s ,C) +s X c(r s ,() +w R (r s ,r c )+e M (r s ), 
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All equations throughout this paper are expressed in Rydberg atomic units. Here t s and e xc are the mean nonintcr- 
acting kinetic energy and the exchange-correlation energy per particle, respectively. For e c we use the Perdew-Wang 
parametrizationJia Here, £m is the average Madelung energy, £m = — 9z/5ro, and ro is the radius of the WS sphere, 
We set for monovalent metals z — 1, and for polyvalent metals we set z* = 1 (for details see Ref.f|). 
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In Eq.|), C = (n T - n x )/(n T + n L ) in which n T and n l are the spin densities of the homogeneous system with total 
density n =F-pt + n l- The quantity wr is the average value (over the WS cell) of the repulsive part of the Ashcroft 
empty coreLLS pseudopotential, 
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and is given by wr — 3r^./r^ where, z is the valence of the atom, 6{x) is the ordinary step function which assumes 
the value of unity for positive arguments, and zero for negative values. 

The core radius is fixed to the bulk value, , by setting the pressure of the unpolarized bulk system equal to zero 
at the observed equilibrium density fi = 3/47Tff : 
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The derivative is taken at fixed r c , and the solution of the above equation gives rf as a function of r s 

_3/2 ( r „ "1 V 2 
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Here, r s = r s (( = 0) is the observed equilibrium density parameter for the unpolarized bulk system, and takes the 
values of 2.07, 2.19, 3.28, 3.99, 4.96, 5.63 for Al, Ga, Li, Na, K, Cs, respectively. Inserting the rf from Eq. (0) into 
Eq. (pi) , the equilibrium r s value of the polarized bulk system is obtained by the solution of the equation 
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= 0. (8) 



The derivative is taken at fixed C & n d rf . The solution gives the equilibrium density parameter f s (£) as a function 
of the polarization. In this procedure, by taking a constant value for the core radius of the pseudopotential, we have 
assumed that the core region of an atom is rigid and does not change in the process of the spin polarization of the 
delocalized valence electrons. This approximation works well when the distance between the neighboring atoms is 
sufficiently larger than the extension of the core electron orbital wave functions. 
The SJM energy for a spin-polarized system with boundary surface isu 



£sjm [n T ,n x ,n + ] = E JM [n T , n h n+] + (e M (r s ) + w R (r s , rf )) / drn + (r) 

+{5v) ws (r s ,r* ) / dr 6(r) [n(r) - n+(r)] , (9) 



+ 1 I dr0([n,n+];r)[n(r)-n + (r)] (10) 



where 

E m [n^,ni,n+] = T s [n T ,n x ] + E xc [ra T , raj 
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2 

and 

^([n,n + ];r)= 2 /^M£W£01. 

Here, n — n-\ + ra x and O(r) takes the value of unity inside the jellium background and zero, outside. The first and 
second terms in the right hand side of Eq.(|To|) are the non-interacting kinetic energy and the exchange-correlation 
energy, and the last term is the Coulomb interaction energy of the system. The effective potential, used in the self- 
consistent KS equations, is obtained by taking the variational derivative of the SJM energy functional with respect 
to the spin densities as 

<ff([^,n l ,n + ];r) = ^—^(E s . m - T.) 

= 4> ([n, n+] ; r) + < c ([n T , nj ; r) + 9(r)(^) ws (r S) rf ), (12) 
where a =|, |. By solving the KS equations 

(V a +< // (r))^f(r)=ef0f(r), a=T,|, (13) 



n(r) = M^, (14) 
ntT (r)= ICWI 2 , (15) 

i(occ) 

and finding the self-consistent values for ef and (f>J, one obtains the total energy. 
In our spherical JM, we have 

Mr) = ^0(R-r) (16) 

in which R = (ziV) 1 / 3 ^ is the radius of the jellium sphere, and n(r) denotes the electron density at point r in space. 
The quantity (<5i>)ws is the average of the difference potential over the Wigner-Seitz cell and the difference potential, 
Sv, is defined as the difference between the pseudopotential of a lattice of ions and the electrostatic potential of the 
jellium positive background. Using the Eq. (21) of Ref. |pj, this average value is given by 
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<^>ws(r a> rf) = -Ui--— . (17) 

Applying Eq. (||) to a metal cluster which contains iVj spin- up, iVj_ spin-down and (= +N^) total electrons, 
the SJM energy becomes a function of N, ( = (JVf —N±)/N, r s , and r,f. The equilibrium density parameter, F a (iV, C), 
for a cluster is the solution of the equation 

= 0. (18) 

r s =f s (JV,C) 

Here again, the derivative is taken at fixed values of AT, £, and rf. For an A-electron cluster, we have solved 
the KS equations^ self-consistently for various spin configurations and r s values and obtained the equilibrium density 
parameter, r s {N, (), and its corresponding energy, E(N, C) = E(N, C, f s (N, () 7 rf) for each allowed spin configuration. 
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-E(N,t,r s ,r*) 



III. RESULTS 

In Fig. [I] we have shown the rf values obtained from the application of Eq. (^) to different r s values. Different 
metals are specified by the rigid squares. It is seen that the plot shows a linear behavior for relatively low electron 
density metals. Inserting these values of rf into Eq. ([[]) and solving the Eq. (||) for different values of polarization, 
C, we have obtained the equilibrium r s values for different bulk metals at various polarizations. The values obtained 
at C = 1 are 2.62, 2.72, 3.69, 4.36, 5.28, and 5.93 for Al, Ga, Li, Na, K, and Cs, respectively. The result for all 
polarizations (0 < ( < 1) is shown in Fig. |^. In this figure, we have plotted the changes in the equilibrium r s values 
relative to the unpolarized value as a function of the polarization for different metals. The results show an increasing 
behavior for all metals but, the increasing rate is higher for higher electron density metals. This property can be 
easily explained in terms of the Fermi holes arround the electrons seen by other electrons with the same quantum 
numbers. In the case of very low electron density metals, the average distance between the electrons is much more 
larger than the effective range of the Fermi hole. Therefore, if some or all of the spin-down electrons undergo a spin 
flip and change to spin-up electrons, the number of spin-up electron Fermi holes will increase. But, these Fermi holes 
will not have any overlap and therefore, effectively no repulsive force will act on the electrons to increase the volume 
of the system. On the other hand, if the average distance between the electrons in a metal be much less than the 
range of a Fermi hole, then these spin flips cause the Fermi holes to overlap and as a result a large repulsive force 
will act on the electrons to increase the volume of the system. This explains the higher increasing rate for Al and 
the lower increasing rate for Cs metals. Figure || shows the plot of (<5v)ws (fs,ff (^s)), obtained from Eq. ( p"7| ) and 
Eq. (^), as a function of the bulk equilibrium WS radius, f s . The rigid squares correspond to different metals. The 
value of (Sv)ws for Na is vanishingly small because, the equilibrium r s value of sodium is very close to 4.18 at which 
this average vanishes. For values of r s greater than 4.18, the correction to the KS effective potantial is positive [see 
Eq.(|l2|)] and this decreases the well depth and cause the electrons relax outward to make the pressure on the system 
(due to jellium model) to vanish. In this case, the leakage of the electrons across the jellium boundary surface is 
increased relative to the simple JM case. On the other hand, for r s < 4.18 the potential depth is increased by the 
correction and the electrons should decrease their relative mean distance to make the pressure on the jellium system 
to vanish. Therefore, the leakage of the electrons decreases relative to the simple jellium model case. In Fig. |^, using 
the Eqs. (0), (|§|), and (|l7] ) we have plotted the variation of (<5v)ws (^s(C)i r f)) as a function of £ for different metallic 
densities. For low electron density metals, it has a decreasing behavior and in the case of K metal, there is a change 
in the sign at high polarizations. However, for Al and Ga, it shows rather different behaviors. That is, all Cs, K, 
Na, Li show a decreasing behavior over the whole range < C < 1 whereas, Al and Ga have decreasing behaviors at 
lower polarizations and increasing behaviors at higher polarizations. These different behaviors can be explained by 
the detail analysis of the Eq. ( |l7|) . The right hand side of Eq. ([I?]) has a physical minimum at r s = y/TErf . In order to 
realize this minimum for a metal over the range < C < 1, the core radius of that metal should satisfy the inequality 

fJQ) R rJl) , , 

<rf < -^=Jr. (19) 
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In the above inequatity, f s (0) and f s (l) are the equilibrium values at £ = and £ = 1, respectively. Examining for 
different metals show that, out of the six mentioned metals, only Al and Ga with respective core radius values of 0.56 
and 0.65 satisfy this constraint. Comparing this figure with Fig. 1 of Ref. |Q, rejzeals the different behaviors for high 
electron density metals, predicted by the stabilized spin-polarized jellium modeltj (SSPJM). The different behavior 
in the SSPJM has roots in the rough estimation used there for the increment in r s due to polarization. But here, we 
could calculate the increment for each metal separately, as shown in Fig. 0. 
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In order to find the equilibrium size of an N electron cluster in the SJM for various spin configurations of Nf— N± = 0, 
iV T - Nt = 2, Nj - Ni = 4, • ••, Nf — Ni = iVf, keeping the total number of electrons, iV T + iV x = N, fixed; we 
solve the Eq. (GJh using the set of self-consistent Eqs. ([l2|)-(^5|). In the above we have assumed an even number of 
total electrons, N. For an odd number of electrons, the differences (JVj — iVj) would also be odd numbers. In Figs. 
^(a)-(c) we have plotted the equilibrium r s values, r s (N, C), obtained for different spin configurations, as a function 
of ( = (jVj — N±)/N for Cs, Na, and Al clusters. To clarify the different behaviors of the closed-shell and open-shell 
clusters, we have studied the N—8, 20, 40 cases which are closed-shell clusters and N—27 which is an open-shell cluster. 
We have also compared the result with the bulk case. In Figs. |j| the dashed lines correspond to the equilibrium WS 
radius of the unpolarized bulk metal, f s (£ = 0). In all the closed-shell clusters, f s (iV, C) is an increasing function of 
C while for the open-shell cluster, N = 27, it is decreasing over < £ < 7/27, and increasing over 7/27 < £ < 1. 
Here, Co = 7/27 corresponds to an electronic spin configuration for which Hund's first rule is satisfied. That is, the 
configuration in which the up-spin shell with / = 3 is half filled. This difference between closed-shell and open-shell 
clusters can be explained as follows. For an open-shell cluster if one increases the spin polarization from the possible 
minimum value consistent with the Pauli exclusion principle, one should make a spin-flip in the last uncomplete shell. 
Because of high degeneracy for the spherical geometry, this spin-flip in the last shell does not change appreciably the 
kinetic energy contribution to the total energy but changes appreciably the exchange-correlation energy which in turn 
gives rise to a deeper effective potential that makes the KS orbitals more localized and therefore a smaller size for the 
cluster. On the other hand, in closed-shell clusters increasing the polarization is always accompanied by transition 
of electrons to unoccupied shells that have larger kinetic energies which leads to larger cluster sizes. In Figs. || we 
notice that as the size of the cluster increases, the plot of f s (N,Q resembles much more to the bulk function f s (C) 
and approaches to it so that, as is expected, we have 

lim f a (N,Q=f s (0- (20) 

N— >oo 

As we see, in all neutral clusters, the inequality 

f s (N,()<f s (Q) (21) 

always holds for the ground state of the cluster in which C = Co- This effect is called self-compressior0 and is due 
to surface tension. Now if we increase the polarization of the cluster, £ = (iVj — N±)/N relative to Co arid obtain the 
equilibrium WS radius by solving the Eq. (|l8|) , we see that beyond some polarization, Ci, the inequality in Eq.(|2~l|) 
changes the direction and the equilibrium r s value of the.-cluster exceeds the bulk value f s (0). This is called self- 
expansion which is also observed in charged metal clusters.E3 This value of Ci is shifted toward zero as the size of the 
cluster increases. Also, comparison of Figs, ^(a)-(c) show that for two clusters with the same N, the value of Ci is 
smaller for higher electron density metal. 

In Figs. ||(a)-(c), we have plotted the quantity (fo)ws for Cs, Na, Al clusters of different sizes, as a function of C, 
using the equilibrium values f s (N, C) in the Eq. ([17|) ; and have compared with their respective bulk functions. It is 
seen that for metals which the inequality in Eq. (EG) does not hold, this quantity has a decreasing behavior whenever 
f s (N, C) has an increasing behavior, and vice versa. For Cs clusters in Fig.||(a), the quantity has positive values over 
the whole range < £ < 1. In Fig. ||(b), for Na clusters, this quantity changes sign beyond some C value. The Al 
clusters in Fig. ^(c) have negative values for (<5v)ws and they show minima as discussed before. 



IV. SUMMARY AND CONCLUSION 



In this work, we have calculated the equilibrium r s values of different metals as a function of their electronic spin 
polarizations, using the stabilized jellium model along with the local spin density approximation. Our calculations 
show an increasing behavior for the bulk Wigner-Seitz radius of electron as a function of polarization. Also we have 
shown that the increasing rate is higher for higher electron density metals. Calculation of the equilibrium r s values for 
closed-shell clusters show similar behaviors as their bulk metals i.e., they also show increasing behaviors as functions 
of C- But, the situation is somewhat different for open-shell clusters. The open-shell clusters show decreasing behavior 
at lower polarizations and increasing behaviors at higher polarizations. The equilibrium r s values of the ground state 
configuration of the clusters are always smaller than the bulk value. This self-compression is due to the surface 
tension. On the other hand, at higher polarizations, the equilibrium r s values exceeds the bulk value r s (0), and this is 
called self-expansion. In conclusion, the SJM, has provided a method which can be used to calculate the sizes of the 
simple metal clusters with minimum possible efforts. More realistic results for open-shell clusters are possible when 
the spherical geometry for the jellium background is replaced by the spheroidal or ellipsoidal shapes. Work in this 
direction is in progress. 
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FIG. 1. The pseudopotential core radius, rj? , in atomic units as a function of the equilibrium r s value of the bulk metal. 
The rigid squares specify different metals. The plot shows a linear behavior at rather low densities. 

FIG. 2. The changes in the equilibrium WS radius relative to the unpolarized value, [f s (C) — *s(0)], in atomic units, for 
different metals as functions of the polarization, £. All have increasing behaviors and the increasing rate is higher for higher 
electron density metals. 

FIG. 3. The average value of the difference potential, in Rydbergs, as a function of the equilibrium WS radius for a bulk 
jellium system. The rigid squares correspond to different metals. 

FIG. 4. The average values of the difference potential, in Rydbergs, as functions of the polarization for different bulk metals. 
The plots show decreasing behaviors for Cs, K, Na, Li over the whole range; whereas the plots for Ga and Al show decreasing 
behaviors at low polarizations, and increasing behaviors at higher polarizations with minima in between. 

FIG. 5. The equilibrium WS radius, f s {N,Q, in atomic units, as functions of the polarization for (a)- Cs, (b)- Na, (c)- Al 
clusters. In all closed-shell clusters (here, A?=8,20,40) f s (N,£) is an increasing function of £ while for the open-shell cluster 
(here, N = 27) it is decreasing over < £ < 7/27, and increasing over 7/27 < £ < 1. The solid lines correspond to the bulk 
f a (£), and the dashed lines correspond to the equilibrium value of the unpolarized bulk system, r s (0). 

FIG. 6. The average values of the difference potential, in Rydbergs, as functions of the spin configurations, Q. In (a)- Cs and 
(b)- Na, for the closed-shell clusters (here, iV=8,20,40) it shows a decreasing behavior over the whole range, the same as the 
bulk case which is specified by a solid line. But, in open-shell clusters (here, N = 27) it is increasing over < £ < 7/27, and 
decreasing over 7/27 < £ < 1. In (c)- Al, the same behaviors hold for lower polarizations but, as discussed in the text, it shows 
a minimum at higher polarizations. 
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